Jamming in Systems With Quenched Disorder 
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We numerically study the effect of adding quenched disorder in the form of randomly placed 
pinning sites on jamming transitions in systems that jam at a well defined point J in the clean 
limit. Quenched disorder decreases the jamming density and introduces a depinning threshold. The 
onset of a finite threshold coincides with point J at the lowest pinning densities, but for higher 
pinning densities there is always a finite threshold even well below jamming. We find that proximity 
to point J strongly affects the transport curves and noise fiuctuations, and observe a change from 
plastic behavior below jamming, where the system is highly heterogeneous, to elastic depinning 
above jamming. Many of the general features we find are related to other systems containing 
quenched disorder, including the peak effect observed in vortex systems. 

PACS numbers: 64.60.Ht,83.60.Bc,83.80.Fg,74.25.Wx 



When a collection of particles such as grains is at low 
densities with little grain-grain contact, the system acts 
hke a liquid in response to an external drive. At higher 
densities, however, significant grain-grain contacts occur 
and the system responds like a rigid solid, exhibiting a 
jamming transition where the grains can become stuck 
under an external drive [H The jamming density is 
termed point J in simple systems such as bidisperse disk 
assemblies [ll-|3|- An already jammed system can be un- 
jammed by the application of shear [5|-|7[ and numerous 
studies have focused on understanding jamming for var- 
ied grain shapes Q, interactions 0, temperatures [lo| . 



and external drives ll|, |12|. Depinning is another ex- 



ample of a transition from a stuck or pinned state to a 
flowing state under an applied drive, and occurs for col- 
lectively interacting particles in quenched disorder such 



as vortices in type-II superconductors |13l - ll5| . colloids in- 



teracting with random or periodic substrates |16l4l8l | , and 
charge-density waves [3]. Above depinning, the particles 
pass from a stationary solid state into either a flowing 
solid or a fluctuating, liquidlike state 1^ IS 16 1. Under- 
standing how quenched disorder affects jamming and how 
jamming-unjamming transitions are related to depinning 
would have a great impact in both fields. 

The first proposed jamming phase diagram for loose 
particle assemblies had three axes: inverse density, load, 
and temperature [1]. Here we propose that quenched 
disorder can form a fourth axis of the jamming phase 
diagram, and show that if a system has a well defined 
jamming transition in the absence of quenched disorder, 
proximity to point J is relevant even for strong quenched 
disorder. We also show that jammed or pinned states 
below point J show profoundly different behaviors in re- 
sponse to an external drive compared to states above 
point J. We find that for varied amounts of disorder, this 
system exhibits many features found in vortex matter 
ll[23,|2l| including a peak effect near point J, suggest- 
ing that jamming may be a useful way to understand 



many of the phenomena found in systems with pinning. 
We study a two-dimensional (2D) bidisperse disk system 
with equal numbers of disks of each size that is known to 
exhibit jamming at a well defined density for zero tem- 
perature and load. The density is defined as the fraction 
(f) of the system area that is covered by the disks. For a 
disk radius ratio of 1 : 1.4, in 2D the onset of jamming 
occurs at « 0.844 [3, i, [III, HI] . Since point J in this 
system is well defined in the absence of quenched disor- 
der, we can determine how jamming changes when we 
add a small amount of quenched disorder in the form of 
randomly placed pinning sites. We also focus on distin- 
guishing the effect of jamming from that of depinning. 
This is particularly important since even non-interacting 
particles that do not exhibit a jamming transition in the 
absence of quenched disorder can still exhibit a finite de- 
pinning threshold in the presence of disorder. 

Simulation- We consider a 2D system of size LxL with 
periodic boundary conditions in the x and y-directions 
containing N disks that interact via a short range repul- 
sive spring force. The sample is a 50:50 mixture of disks 
with interaction ranges rA and tb, where ta = 1.4r's. 
Nj is the number of particles in the sample at the jam- 
ming transition. To initialize the system, we place disks 
in nonoverlapping positions and then shrink all disks, add 
a few additional disks, and reexpand all disks while ther- 
mally agitating the disks until reaching the desired den- 
sity. We employ overdamped dynamics where the equa- 
tion of motion for the disks is rjdRi/dt — X^j^itj ^(^cff ~ 
|rij|)+Fp+F£i. Here the damping 



J, where r^Q) 



|r,,|)(r,,/|r,,|)e(i?^ 

constant r/ ^ 1, k = 2uu, and J^^ff = fi + i" 
is the radius of disk i{j). The driving force F^i = -Fbx is 
applied to all the disks uniformly. The pinning force 
is modeled as arising from Np non-overlapping attractive 
parabolic traps with a maximum force of Fp and a cut- 
off radius of r^/2 to ensure that only one grain can be 
trapped by any one pinning site. In the absence of other 
grains, an isolated grain will depin when Fo > Fp. In 



2 



•z. 

> 

V 

1 



(a) ' 1 


1 / 


(|)/(|) =0.761/ 
J\ 1 


(|)/(|)=0.014- 




0.02 0.04 0.06 0.08 0.1 




FIG. 1: (Color online) (a,b) The normalized average grain 
velocity {Vx)/N vs external drive Fd in samples with Fp — 
2.0. (a) Np/Nj = 0.415. Right (red) curve: at (l)/(j>.j = 0.014 
depinning occurs at Fc — 2.0. Left (black) curve: at (p/cpj — 
0.761, Fc is much lower, (b) Np/Nj = 0.09267. Lower (black) 
curve: at </>/</>j = 0.947, there is no finite depinning threshold. 
Upper (red) curve: at (t>/(t>j = 0.99, Fc is finite, (c) Disk 
positions in the pinned state for the system in (a) with — 
0.761. (d) The disk trajectories over a period of time for the 
system in (c) at Fd ~ 0.25 showing plastic flow above the 
depinning threshold, (e) Disk positions in the pinned state 
for the system in (c) with = 1.03878. (f) The disk 

trajectories over a period of time for the system in (e) at 
Fd = 1-1 Fc showing elastic flow. 



the absence of pinning, this model has been well studied 

and has a jamming density of 4>j ~ 0.844 @, im. To 

determine if the system is pinned, we measure the total 
average grain velocity (V^) — X^i^i • x as a function of 
Fo, where is the velocity of grain i. For most of the 
results presented, Fo is slowly increased in increments 
of SFd = 5 X 10~® and we wait 5 x 10^ simulation time 
steps after each force increment to ensure that the system 
reaches a steady state response. 

In Fig. l(a,b) we plot {Vx)/N versus Fd for differ- 
ent pinning densities and particle densities, measured in 



terms of Nj and i^j. At Np/Nj = 0.415 and (p/(pj 
0.014 in Fig. 1(a), the system is pinned for drives up 
to Fd = 2.0. At this low density, there are no inter- 
actions between the disks and the depinning threshold 
Fc is solely determined by the pinning force. For the 
same pinning density at = 0.761, there are more 

disks than pinning sites; however, as shown in Fig. 1(a), 
there is still a finite depinning threshold of Fc = 0.27. 
Since not all of the disks can be captured by pinning 
sites, the existence of a finite depinning threshold indi- 
cates that the disks trapped by pins must be blocking the 
flow of the disks that are not in pins. Thus, in pinned 
states such as that illustrated in Fig. 1(c), some form 
of jamming must be occurring. The depinning thresh- 
old can vanish when the pinning density is reduced, as 
shown in Fig. 1(b) for a sample with Np/Nj = 0.09267 
and = 0.947. Here Fc = and the velocity re- 

sponse is nonlinear. The depinning threshold remains 
zero as (j)/4>j is reduced until N drops below Np and 
all the disks can be trapped, giving Fc = Fp. At high 
disk densities such as 0/0,7 = 0.99 in Fig. 1(b), Fc is 
finite. Figure 1(b) also shows an interesting crossing of 
the velocity-force curves, caused by the sudden jump in 
{Vx)/N at depinning for the = 0.99 sample. The 

depinning is elastic for 4>/(f>j = 0.99, so all the disks begin 
to move simultaneously at the same velocity with only 
small localized rearrangements of the disk packing. In 
contrast, for (f>/4>j = 0.947 the depinning is plastic. Some 
disks remain pinned while other disks move past, and the 
pinned disks do not depin until a much higher Fd is ap- 
plied (not shown) . An example of the plastic flow occur- 
ring at depinning appears in Fig. 1(d) for a sample with 
Np/Nj = 0.415 and 0/0j = 0.761 at Fd = 0.25, above 
the depinning threshold. Certain disks are always pinned 
while rivers of disks flow around them. In the jammed 
state, the disk packing acts like a solid and a small num- 
ber of pinning sites can trap all the disks, as in the sam- 
ple with Np/Nj = 0.09267 and (/-/0j = 0.99 in Fig. 1(b). 
The depinning from the jammed state is elastic, as illus- 
trated for (p/ipj = 1-03878 in Fig. l(e,f). Below jamming 
in samples containing a small number of pinning sites, 
such as the Np/Nj = 0.09267 and 0/0j = 0.947 sample 
shown in Fig. 1(b), some but not all of the disks can be 
immobilized, so the depinning is plastic. 

Figure 2(a) shows Fc versus 4>/(j)j for a series of dif- 
ferent values of Np/Nj, and Fig. 2(b) shows a blowup of 
the same data near 0/0j = 1.0. For Np/Nj > 0.138 the 
depinning threshold is finite for all (t>/4>j- For Np/Nj < 
0.138 at low 0/(/)j, the depinning force reaches its maxi- 
mum possible value of Fc = Fp = 2.0 since all the disks 
can be trapped by pinning sites without any overlap. 
There is a region of intermediate 4>/(f>j where there is 
no depinning threshold. At higher 4>/(j)j, Fc becomes fi- 
nite again and the system undergoes plastic depinning. 
The value of Fc increases with increasing 4>/(j)j in this 
regime until Fc peaks near 4>/4>j = 1.0 when the system 



3 




FIG. 2: (Color online) vs j for Np/Nj = 0.828 (Q), 
0.415 (red ■), 0.277 (green 0), 0.138 (blue A), 0.09267 (T), 
0.0346 (red >), 0.00692 (green +), and 0.00138 (blue x). (a) 
A log-linear plot over the full range of (?!)/(/) j shows that the 
depinning threshold is always finite for Np/Nj > 0.138, while 
for Np/Nj < 0.138 there is a pinned regime at low 4'/4'J when 
all the disks are trapped by pinning sites and another regime 
of finite Fc at high 4'/4'J where jamming occurs, (b) A blowup 
of the region near (f>/<pj = 1.0 shows that the onset of jam- 
ming is marked by either a drop in F^ for the higher pinning 
densities or a peak in Fc for the lower pinning densities. 
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FIG. 3: (Color online) Black circles: The value of Np/N,j at 
which a finite depinning force first appears vs <j). Red squares: 
The value of Np/Nj at which Fc reaches its peak value vs (p. 
The system depins elastically for <jf> just above the peak value 
of Fc- For small <j), all of the particles are directly pinned at 
pinning sites, while for close to = 0.844, the particles 
are pinned due to jamming. 

jams. At jamming, in principle a single pinning site could 
pin the entire granular packing; however, we find that Fc 
at jamming decreases as the pinning density is reduced. 
At and above jamming, the depinning becomes elastic 
and Fc drops with increasing (j)/(t>j due to the increasing 
stiffness of the jammed solid. 

Using velocity force curves, we identify the onset of a 
finite depinning threshold as a function of (j) for different 
values of Np/Nj, as well as the value of cj) at which Fc 
reaches its peak value for (p near (pj . The result is shown 
in Fig. 3. There is a dome feature connecting the low 



(j) limit, where the disks are directly pinned at pinning 
sites, with the high limit, where the pinning occurs 
when the grains become jammed. Near we find that 
the quenched disorder density Np/Nj can be considered 
as a new axis of the jamming phase diagram, and that 
the jamming density decreases with increasing pinning 
density. Although the onset of jamming can be defined 
as coinciding with the appearance of a finite depinning 
threshold, it could also be defined as coinciding with the 
transition from plastic to elastic depinning, where Fc 
reaches its peak value. Figure 3 shows that these two 
definitions are not identical but track each other closely 
in the region near . The onset of elastic depinning con- 
tinues to produce a peak in Fc up to the highest pinning 
densities we have considered. 

We can make a simple argument for how quenched dis- 
order reduces the jamming density. The average distance 
between pinning sites is Ip cx pp where pp = Np/L'^. 
We can estimate a correlation length ^ by assuming 
that ^ grows as jamming is approached according to 
^ cx: {(j>j — (j))^'^ ■ Jamming should occur when Ip = ^, 
or when pp oc {(f) j — (f)f''^ ■ A fit of the onset of finite Fc 
for (f) > 0.8 or a fit to the peak value of Fc both give 
a linear dependence on 0, implying that v = 0.5. This 
value is close to some threshold predictions [11]; how- 
ever, caution must be taken in comparing our exponent 
to systems without quenched disorder, since the pres- 
ence of quenched disorder or the fact that we are driving 
our system could fundamentally change the nature of the 
jamming compared to the disorder free case, and addi- 
tional corrections to scaling could be relevant [^j j. 

The onset of jamming can also be detected by an- 
alyzing the velocity noise fluctuations using the power 
spectrum S{f) obtained from the time series of the av- 
erage disk velocities, S{f) = \J exp(-i27r/t)(14)(t)dtp. 
The noise power 5*0 is defined to be the integrated 
noise power in the lowest octave of the spectrum. We 
take Fd = l.li^c- Figure 4(a) shows Sq versus (j) for 
Np/Nj = 0.277 and Fig. 4(b) shows the same quantity 
for Np/Nj = 0.0346. For Np/Nj = 0.277, the noise 
power So is low for < 0.3 since there are very few 

collective interactions at these low densities that could 
give rise to low frequency noise. Sq decreases rapidly 
with increasing (p/cpj for (f)/4>j > 0.95 after the system 
jams and transitions to elastic depinning. At interme- 
diate (p/cpj, large velocity fluctuations occur and pro- 
duce a broad band noise signal, as shown in Fig. 4(c) 
for (t)/(t>j = 0.968 where the depinning is plastic. The 
solid line in Fig. 4(c) is a flt to 1//°'^. The appearance 
of 1// type noise is known to be associated with plastic 
depinning [l^[23|- For (p/cj),] — 1.141 in the elastic de- 
pinning regime. Fig. 4(c) shows that the noise power is 
considerably reduced and the spectrum has a peak at fl- 
nite frequencies with several higher harmonics, indicative 
of a narrow band noise signal. The appearance of narrow 
band noise in driven systems with quenched disorder is 
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FIG. 4: (Color online) (a,b) The noise power So obtained 
from the time series of the velocity fluctuations vs (f>/(j>j for 
Fd = 1.1-Pc. (a) Ny/Nj = 0.277. So drops at the onset of 
jamming and at low 4>/4>J- (b) Np/Nj — 0.0346. There is a 
pronounced peak in just below jamming, (c) The power 
spectrum S{f) vs / from the system in (a) for (upper red 
curve) <p/<pj ~ 0.968 where the depinning is plastic and (lower 
blue curve) 4'/4'J ~ 1-141 where the depinning is elastic and 
a narrow band noise signal appears. The solid black line is a 
flt to 1//"'^. (d) S{f) vs / for the system in (b) for (upper 
red curve) 0/(/>j — 0.989 and (lower blue curve) (j>/(j>j = 1.141 
showing narrow band noise in the jammed phase. The solid 
black line is a fit to 1//^. 



associated vi^ith the formation of a moving solid [l9j, |25 1 
and is commonly referred to as a washboard frequency. 
This is consistent with the moving jammed packing act- 
ing like a rigid solid. For the smaller pinning density of 
Np/Nj = 0.0346, Fig. 4(b) shows a pronounced peak in 
So at = 0.989, which corresponds to a density just 

below the peak in Fc. Sq then drops off rapidly as the sys- 
tem enters the jammed phase and depins elastically. In 
Fig. 4(d) we plot the power spectra for 4>/(l)j = 0.989 and 
1T41 for Np/Nj = 0.0346. There is broad band 
noise in the plastic flow regime at (p/tpj — 0.989 and nar- 
row band noise in the elastic flow regime at (f>/(j>j — 1.141. 

The depinning-jamming transition we observe has 
many similarities to the peak effect phenomenon in type- 
II superconductors with vortices moving through random 
disorder. The peak effect can occur as a function of vor- 
tex density. At low densities, the depinning threshold is 
high since vortices can be pinned individually. At inter- 
mediate densities, the depinning threshold remains at a 
low value until some higher density is reached where there 
is a significant increase in the depinning threshold to a 
peak value. This is associated with a large increase in the 
noise power 13, 2^ 21 1. The peak effect and the noise 
features become more prominent for cleaner samples with 
less pinning [2l|. All these features are captured in our 
results. The standard interpretation of the peak effect is 
that it marks a transition from a weakly pinned solid to a 



more strongly pinned disordered state. Our results sug- 
gest that the peak effect may be a general phenomenon 
in other systems with quenched disorder close to some 
type of phase transition. 

In summary, we show how jamming behavior changes 
with the addition of quenched disorder using a simple 
model of bidisperse disks that exhibit a well defined jam- 
ming density in the absence of quenched disorder. 
We propose that quenched disorder represents a new 
axis of the jamming phase diagram and that increasing 
quenched disorder density decreases the disk density at 
which the system jams. At low disorder densities, the 
disk density at which a finite depinning threshold ap- 
pears coincides with point J. There is also a reentrant 
finite depinning threshold at low disk densities when all 
the disks are directly pinned. We find a maximum in 
the depinning threshold just at the onset of jamming for 
low disorder densities. When the disorder density is suffi- 
ciently large, the depinning threshold is finite for all disk 
density values; however, proximity to produces clear 
effects in the form of features in the velocity force curves 
as well as noise fluctuation signatures. Below jamming, 
the depinning is characterized by plastic flow and 1/f 
noise characteristics with strong heterogeneities in the 
system, while above jamming, the depinning is elastic 
with all the particles moving together and is character- 
ized by a washboard noise. For high disorder density the 
onset of jamming is associated with a drop in the depin- 
ning threshold as opposed to the peak in depinning found 
at low disorder density. Our results show many similar- 
ities to the peak effect observed in high-temperature su- 
perconductors where a peak in the depinning threshold 
occurs at both low and high vortex densities. Our re- 
sults should be relevant for systems exhibiting depinning 
transitions and jamming. 
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